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ABSTRACT 

We consider the consequences of appreciable line optical depth for the profile shape of X-ray 
emission lines formed in stellar winds. The hot gas is thought to arise in distributed wind shocks, 
and the line formation is predominantly via collisional excitation followed by radiative decay. 
Such lines are often modelled as optically thin, but the theory has difficulty matching resolved 
X-ray line profiles. We suggest that for strong lines of abundant metals, newly created photons 
may undergo resonance scattering, modifying the emergent profile. Using Sobolev theory in 
a spherically symmetric wind, we show that thick-line resonance scattering leads to emission 
profiles that still have blueshifted centroids like the thin lines, but which are considerably less 
asymmetric in appearance. We focus on winds in the constant-expansion domain, and derive an 
analytic form for the profile shape in the limit of large line and photoabsorptive optical depths. 
In this limit the emission profile reduces to a universal shape and has a centroid shift of —0.24^00 
with a HWHM of 0.63woo. Using published data for Chandra observations of five emission lines 
from the O star £ Pup, we find that the observed HWHMs are somewhat smaller than predicted 
by our theory; however, the centroid shifts of all five lines are consistent with our theoretical 
result. These optical-depth effects can potentially explain the more nearly symmetric emission 
lines observed in £ Ori, 9 1 Ori C, and 5 Ori by Chandra, although an alternative explanation is 
required to account for the unshifted peak line emission. We also consider enhanced reabsorption 
by continuous opacity as line photons multiply scatter within an optically thick line, and find for 
lines with optical depths of a few that such reabsorption can further reduce the line asymmetry. 
It also reduces the line equivalent width, but probably not enough to alleviate the problem of 
sub-solar metallicities inferred from O star X-ray spectra by ASCA, unless the width of the 
resonance regions are superthermally enhanced. 
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1. INTRODUCTION 

Observations by the Chandra and XMM-Newton 
telescopes are providing first-ever resolved emission 
line profiles from early-type stellar winds in the 
X-ray band (the O star C Pup by Kahn et al. 2000 
and Cassinelli et al. 200f ; the O star 9 1 Ori C by 
Schulz et al. 2001; the O star C Ori by Waldron 
& Cassinelli 2001; and the O star 5 Ori by Miller 
et al. 2001). Already, many early- type stars were 
known to be X-ray sources from observations by 
EINSTEIN and ROSAT (Harndcn et al. 1979; 
Seward et al. 1979; Seward & Chlcbowski 1982; 
Pollock 1987; Pollock, Haberl, & Corcoron 1995; 
Kudritzki et al. 1996; Berghofer et al. 1997). 
Although early attempts sought to explain the X-ray 
emission via stellar coronae (e.g., Cassinelli & Olson 
1979; Waldron 1984), the currently favored model is 
based on considerations of stellar- wind shocks (e.g., 
Lucy & White 1980; Lucy 1982; Owocki, Castor, & 
Rybicki 1988; Fcldmcicr, Puis, & Pauldrach 1997). 
The highly successful line-driven wind theory for 
explaining hot star winds (Lucy & Solomon 1970; 
Castor, Abbott, & Klein 1975; Friend & Abbott 
1986; Puis, Pauldrach, & Kudritzki 1986) is well 
suited to generating instabilities leading to shock 
structures distributed throughout the wind flow and 
giving rise to the production of hot gas. 

The confrontation of theory with data has led to 
mixed results. For example, Feldmeier et al. (1997) 
conducted extensive modelling of wind shock X-ray 
production using time-dependent hydrodynamical 
simulations. Their sophisticated treatment includes 
a voluminous line list, but is limited to spherical 
symmetry, and only models for O stars have been 
investigated. Within this context, their efforts reveal 
that a single high-density shell will tend to dominate 
the observed X-ray emission. Such a feature should 
lead to significant variability which has not been 
observed (e.g., see Berghofer & Schmitt 1995). 
Phenomenologically, perhaps their results apply in 
a "piecewise spherical" sense, whereby these X-ray 
bright events no longer occur in a single global shell, 
but as numerous smaller events occuring at different 
times. In this way, the variability can be suppressed 
(e.g., Oskinova et al. 2001). 

Aside from the issues of global X-ray production 
and variability, the resolved emission line profiles in 



four O stars (( Pup, £ Ori, 9 1 Ori C, and S Ori) 
appear to be truly problematic for the current theory. 
The X-ray emission is expected to be dominated 
by a spectrum of optically thin lines. In addition 
to the intervening ISM, the stellar wind itself can 
significantly contribute to the attenuation of shock 
generated X-rays. Consequently, the observed 
emission profiles are expected to be asymmetric, 
with peak line emission appearing blueshifted of 
line center, with the redshifted emission being more 
suppressed because there is a larger column depth 
of attenuation to the receding hemisphere of the 
wind (MacFarlane et al. 1991; Ignace 2001; Owocki 
& Cohen 2001). However, the observed line profiles 
for ( Ori, 9 1 Ori C, and 5 Ori appear to be quite 
symmetric. The emission lines for £ Pup are far more 
consistent with the expectations of distributed wind 
shocks, but even for this star, the line morphologies 
may not be consistent with the current theory. So 
either the underlying wind model, or simply the 
assumptions used in deriving the emission-line profile 
shapes, must be modified. It is this latter case that 
is the focus of this work, with particular attention 
given to £ Pup since it appears to best exemplify 
current standard wind-shock expectations. 

Raymond & Smith (1977) made a major 
contribution to interpreting the X-ray emission from 
hot astrophysical plasmas by providing a cooling 
function that, among other processes, accounts for 
the emission by optically thin lines. Their results, 
along with better line data, have been used in the 
modelling efforts of numerous studies, including 
the topic at hand of stellar winds. However, some 
have questioned the assumption of effectively thin 
emission among the stronger wind lines in the X-ray 
band (e.g., Hillier et al. 1993; Kitamoto & Mukai 
1996; Kitamoto et al. 2000). As will be discussed 
in §2, we find that it is possible for some lines to 
be optically thick. In §3, the formation of optically 
thick line profiles in distributed X-ray emitting 
zones for a stellar wind is considered using standard 
Sobolev theory. In particular, focus is given to the 
case of constant spherical expansion appropriate for 
relatively dense winds in which the X-ray emission 
escapes only from several radii out in the flow. This 
case is germane to the O stars which have been 
observed with Chandra and XMM. Optically thick 
emission not only alters the profile diagnostics, it 
presents the possibility for X-ray photon destruction, 
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which would invalidate the Raymond & Smith 
assumption for those lines. In §4, we consider 
the role of photon absorption in X-ray bound-free 
continua within the resonance zones in which the 
photons are created. A discussion of our results in 
relation to published observations is given in §5. 



2. THE OPTICAL DEPTH OF X-RAY 
LINES 

Should one expect wind emission lines produced 
in the X-ray band to be optically thick? As 
already stated, the primary process by which X-ray 
line photons are created is collisional excitation 
followed by radiative decay. For minority species, 
radiative recombination into excited states may also 
contribute, but for our purposes this may be viewed 
as essentially a collisional process as well. The key 
assumption is that the photon creation rate per 
volume scales with the square of the density. To 
estimate the optical depth, we consider a slab of 
X-ray emitting shocked gas. For a slab of initial 
temperature Tx, constant electron number density 
n e , and volume cooling rate A(Tx), the cooling time 
will be given roughly by 



5 kAT 
2 n c A{T x )' 



(1) 



The AT is to represent the range of temperatures 
over which a particular ion state i of some atomic 
species j exists in the slab for the purpose of 
estimating the line optical depth. The expression 
thus describes the cooling time for a particular ion 
species from T x + 0.5AT to T x - 0.5AT. 



The line optical depth is given roughly as 



r w ni(Til c Ap 1 , 



(2) 



where n\ is the population number density of the 
ground state for the transition of interest, cr; is the 
line integrated optical depth as given by 



m c c 



fl2, 



(3) 



which applies for a two-level atom that is strongly 
NLTE, l c is the cooling length as given by 



lc = Vt c 



(4) 



for V the speed of the postshock material, and 
Av = vd/X is the frequency width of the line opacity 
for Doppler broadening vd and line wavelength A. 

For simplicity the Doppler broadening is just 
taken as the thermal speed of the hot gas, although 
turbulent motions could be important and somewhat 
larger. For NLTE, the population number density 
of the ground state is roughly equal to the entire 
number density for that ion n\. Taking this to be 
the case, ni can be re-expressed in terms of the 
ionization fraction g; of the ion to all ions of the 
species and the relative abundance Aj of the species 
to all nuclei tin, giving n\ « %A)n^. Thus the 
optical depth of the line formed in just this one slab 
becomes 



5 k AT q\ Aj n^aiXV 
2n c A(T x )v th 



(5) 



We estimate the optical depth for O VIII Lya 
at 19.0 A as an example of a strong line in a gas 
with a solar composition of oxygen (Aj = 9 x 10~ 4 ). 
From Sutherland & Dopita (1993), the ionization 
fraction of O VIII peaks at around 2 million K with 
a value of qi « 0.5, which predominantly exists over 
a temperature range of AT w 2 million K. The 
volume cooling rate is about 6 x 10~ 23 erg s _1 cm 3 . 
The necessary shock jump is 500 km s _1 , and the 
thermal speed for hot oxygen is 45 km s _1 . Hence 
the optical depth for this one slab becomes t « 10 2 . 
This demonstrates that it is entirely reasonable that 
the strong lines of metals with Aj ;> 10~ 4 , which 
includes many of the lines observed in hot star X-ray 
spectra, could be optically thick. 

Having shown that optical thickness may be 
relevant for strong resonance lines, what should 
we expect for its influence on the emission profile 
shape? The observed profiles for £ Pup (Kahn et al. 
2001; Cassinelli et al. 2001) are asymmetric with 
blucshiftcd peak emission. As will be discussed 
further at a later point, photoabsorption by the 
ambient wind component generally accounts for 
this asymmetry. However, the effect of line optical 
depth is to alter the primary direction of photon 
escape. For example, the emissivity for a line that is 
thin is isotropic, and so there is no preferred escape 
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direction for photons. For a radially narrow but 
optically thick slab as considered above, photons 
escape primarily normal to the slab face, since the 
optical depth in that direction is least. For a large 
number of slabs at a fixed radius and uniformly 
distributed in solid angle, and ignoring the wind 
photoabsorption, those slabs that are fore and aft 
will tend to be the brightest based on this argument. 
Now accounting for the asymmetric photoabsorption 
across the line frequencies, this suggests that peak 
blueshifted emission should occur near the maximum 
Dopplcr shift. Instead, the observed emission profiles 
tend to have peak emission intermediate between 
line center and the extreme wing. 

One other possibility is that photons escape 
preferentially in the lateral direction. The motivation 
here is that in outflows, velocity gradients can lead 
to just such escape avenues. In this case, consider 
numerous discrete slabs distributed across a spherical 
shell. All else being equal, the brightest slabs will 
be those located near the plane of the sky owing to 
the lateral escape of radiation, and corresponding 
to frequencies around line center, which is more in 
accord with observations. 

However, we caution that the Sobolev 
approximation applies when the Sobolev length Z so b 
is significantly smaller than the cooling length l c 
(i.e., the extent of the hot gas region). A rough 
expression for a ratio of the two is 

where n c is in cm -3 and i> s h is the speed of the 
shocked gas in the observer's frame. We will 
consider a spherically symmetric expanding wind 
with ambient number density 

"wW = -: — TT' ( 7 ) 

47T (i ran r z v T {r) 

where M is the mass-loss rate, v T is the radial wind 
speed, and jj, is the mean molecular weight per 
particle. The electron number density will naturally 
scale with the wind density. For typical hot stars, 
n w ranges from 10 10 -10 12 cm -3 at the wind base 
and drops to 10 7 — 10 9 by about 3i?*, thereafter 
decreasing as r~ 2 . The Sobolev approximation 



should thus be valid at most radii, except possibly 
in the densest regions of the inner wind acceleration 
where v r is relatively small. But even at small 
radii, the Sobolev approach may be adequate if the 
shocked gas is of low density, as would be the case if 
the emission arises in reverse shocks. Alternatively, 
if i?* 3> / S ob l c , then there exist numerous 
pockets of hot gas within a Sobolev length. These 
would likely be optically thin, so treating them as 
statistically homogeneous allows us to recover the 
Sobolev approximation. 

However, to say that the photons escape 
laterally implies that the lateral optical depth is 
less than the radial optical depth. So the effect 
of the velocity gradient carves the slab up into 
numerous radiatively independent and elongated 
"columns" . The short length of these columns is 
the Sobolev length, and the long length is the slab 
thickness. Given that the optical depth through 
the slab's thickness was shown to be w 10 2 , the 
columns can still be thick in the lateral direction for 
Ic/hob ranging from 10-100. This range is entirely 
reasonable given the preceding estimates, and so we 
next apply the Sobolev theory to deriving emission 
profile shapes to model the effects of lateral photon 
escape. 

3. OPTICALLY THICK X-RAY 
EMISSION PROFILE SHAPES 

3.1. Formalism for the Line Formation 

For an envelope in bulk motion such that the 
flow speed greatly exceeds the thermal broadening, 
the locus of points contributing to the emission 
at any particular frequency in the line-profile is 
confined to an "isovelocity zone" (c.f., Mihalas 1978). 
These zones are determined by the Doppler-shift 
formula, namely 

= vo (l - ^) , (8) 

where v 2 is the Dopplcr-shiftcd frequency and 
v z = — v(r) • z is the projection of the flow speed onto 
the line of sight along z. The familiar Sobolev theory 
for line-profile calculation in moving media reduces 
the radiation transfer in the envelope to radiation 



- 5 - 



transfer within isovelocity zones (i.e., the transfer 
between different points in the envelope decouples 
if their relative velocity shift exceeds the local line 
broadening). Thus, the emergent line intensity along 
a given ray through the extended stellar envelope is 
given by 



I v = S v {r) 



1 - e 



-ri(r,n) 



(9) 



where S v is the source function, n is the Sobolev 
optical depth of the line where the ray intersects the 
isovelocity zone, and r c is the intervening absorptive 
optical depth. 

Allowing for line emission from resonance line 
scattering and collisional de-excitation, and defining 
e as the ratio of the collisional de-excitation rate 
(related to the rate of collisional excitation) to that 
of spontaneous decay, the source function is given by 



S v = 



/3 C Iv,* + e B v 
/? + £(!-/?)■ 



(10) 



The two parameters f3 c and (3 are respectively the 
penetration and escape probabilities, defined as 



i 



n 



■da 



(ii) 



and 



4vr J to 



■ da 



n 



(12) 



Adopting standard cylindrical coordinates 
(p, a, z), with the observer at +oo, the line optical 
depth is given by 



for p = cos 9 the direction cosine from the radial to 
the observer. 

For the case of X-ray emission, we assume that 
the star is X-ray faint (i.e., no basal corona) such 
that = at X-ray energies. Also, the destruction 
probability e depends on density, and so drops with 
radius. To make this explicit, we parameterize 



£o p/po, 



(15) 



where po and eo are fiducial values with eo <C 1 being 
typical. Under these conditions, the source function 
becomes 



(16) 



To compute the total emergent flux at a given 
Dopplcr shift in the line profile, we integrate over all 
intensity beams to obtain 



1 

~D 2 



I v {p, z)pdpda, 



(17) 



where D is the distance from the Earth. Self- 
absorption by the hot plasma is ignored. The 
attenuation is therefore entirely from the cool wind 
component intervening between the observer and the 
point of emission. (For actual observations, one must 
also correct for the energy-dependent interstellar 
extinction, but we ignore its effects here since it does 
not alter the profile shape over the narrow extent of 
the line width.) 



The wind optical depth is given by 



k c p w (p, z) dz 



(18) 



n(p,z) 



a i n\ 5(v — v z ) dz' = 



<7i n k (r) A 
\dv z /dz\ p 



(13) 



where cr; is the frequency-integrated line cross- 
section, and the denominator on the right-hand side 
is the line-of-sight velocity gradient. In the case of a 
spherical flow, the line-of-sight velocity gradient is 



dVz = 2^ + / 1 _ 2\V 

dz dr r ' 



(14) 



with opacity k c and density p w . The dominant 
opacity at the X-ray energies is photo-absorption by 
K-shell electrons. This opacity is calculated from 



Kc(E) = Y,—°i{E)- (19) 

p N m H n N 

The opacity is a summation over cross-sections 
presented by different atomic species j and weighted 
by the relative abundance n-Jn^, for tin the number 
density of nuclei. 
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The line source function and a prescription 
for the wind attenuation are all the required pieces 
necessary to compute X-ray emission line profiles 
from stellar winds. The above expressions are 
applied in the next section for the specific case of a 
constant expansion flow. 



that 

/•OO 

t c « 1 = / K c p w dr. (20) 

Jri 

Assuming that v(r) = allows for an analytic 
evaluation of this integral, yielding 



3.2. The Constant Expansion Case 

3.2.1. Justification of Assumption 

The constant expansion wind is the most 
relevant case to consider for the problem at hand, 
since the continuum absorption by the winds of 
most O stars are expected to be optically thick 
over much of the X-ray band. The implication is 
that the observed X-rays emerge from the wind at 
depths of a few stellar radii or more. For the O star 
winds, the standard /3-law velocity prescription for 
the radial distribution of wind speeds begins to 
asymptote around 1 to 1.5 J2* above the wind base 
for the commonly used value of [3 1. Consequently 
for many lines, it is sufficient to model the emission 
profiles as if the entire wind were undergoing 
constant expansion. This case is also convenient 
because the expressions for the radiative transfer 
become analytic and provide interesting insight for 
better understanding the resulting profile shapes. 
However, there is at least one observational example 
of a line that seems to originate from deeper levels 
in the flow where the constant expansion assumption 
no longer applies (e.g., He- like S XV observed in 
C Pup by Cassinelli et al. 2001). In §§, we comment 
on the effect of the wind acceleration zone on the 
profile shape. 

We begin by motivating the constant expansion 
assumption on more quantitative grounds. Owocki & 
Cohen (1999; for OB stars) and Ignace & Oskinova 
(1999; for WR stars) presented a scaling analysis 
for the X-ray emission from hot star winds based 
on an "exospheric" approximation. The observed 
X-ray emission (assumed to arise from an optically 
thin hot plasma) emerges only from radii exterior 
to the optical depth unity surface of radius r\ , with 
X-rays at smaller radii assumed to be completely 
attenuated. The extent of this convenient scale r\ 
derives from the photoabsorption by the wind along 
the line-of-sight. It is determined by the condition 



n(E) = — Kc (E). (21) 

This may seem circular, since we apply the constant 
expansion assumption here to justify that constant 
expansion should be adequate for the emission line 
calculation. But in fact, the integral of equation ( |20| ) 
minimizes the value of r\ in the case of v T = Vqq, 
hence this is the conservative approach to showing 
that r\ ;> 2i?* is valid over much of the X-ray band. 

Using typical O star and wind parameters as 
in §2, a plot of r\ versus wavelength is shown in 
Figure |l|. The "jagged" features correspond to 
prominent photoabsorptive edges. Observed strong 
lines of H and He-like ions in stellar winds are 
generally located between 5 and 30 A, as indicated 
at top. Note that 1 keV corresponds to about 12.4 A 
(indicated by the vertical line in the plot), and so 
harder X-ray energies are at shorter A, and vice 
versa. A horizontal line is provided where r\ = 2i?*, 
beyond which the wind is essentially in constant 
expansion for O star winds. Many of the strong lines 
of interest occur at wavelengths for which r\ ;> 2i?» , 
where the approximation that v r « Vqq will hold. A 
similar figure is provided by Cassinelli et al. (2001) 
in their discussion of Chandra observations of £ Pup. 
They even tabulate r\ values for several lines in their 
Table 2. It happens that ( Pup has a slightly higher 
M and lower than we have used, so that even for 
Mg XII at 8.42 A, Cassinelli et al. find n = 2.8i?*. 

3.2.2. Applying the Constant Expansion 
Assumption 

Concluding that constant expansion is a 
reasonable assumption, we now consider again the 
line formation process under this condition. /From 
equation (|l4|), the line-of-sight velocity gradient 
reduces to 
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n 



<Ji nx(r) Ap 
( Voo /r) (1 - v?) 

r 



T 



(22) 
(23) 



where an optical depth scale T has been introduced. 
This scale is the "integrated optical depth" of the 
line along the line-of-sight (Groenewegen, Lamers, & 
Pauldrach 1989). The parameter T can be related to 
the column density of scatterers N\ as follows, 



m dr 



n 



dv T . 



(24) 
(25) 



Defining T = J„ tj dw for w = w r /woo, we obtain 



(26) 



Unfortunately in the constant expansion 
approximation, the Sobolev optical depth of 
equation (^) contains a pole at fi = ±1, 
corresponding to the extreme blue and red shifts. 
Physically of course, there is no real singularity 
problem. Because the wind is in constant expansion, 
there is no radial velocity gradient, hence the 
actual optical depth in the vicinity of /z = ±1 is 
just (voa/vth)T, for vth the thermal broadening, 
because the radiative transfer corresponds to the 
case of a slab moving at constant speed Voq. Since 
Voa/vth > 10 can be expected, values of T exceeding 
a few will ensure that 77 for ji f» ±1 will be large 
enough that the exponential factors in the Sobolev 
intensity of equation (^) and the escape probability 
of equation ( |l2|) are effectively zero, hence the "pole" 
problem will not be of significant concern. There 
are two subsequent implications for the constant 
expansion wind, (a) Instead of the line-of-sight 
integrated optical depth, T actually corresponds to 
T{ at fi = and is thus a lateral optical depth scale. 
This is reasonable since for constant expansion, 
the line-of-sight velocity gradient arises only from 
the geometrical divergence of the spherical wind 
outflow. This scale is also the minimal Sobolev 
optical depth of 77 (/z) for any fixed r. (b) The 



radial optical depth can be large, and so there is 
the possibility for effective photon trapping and 
substantial reabsorption by continuum opacity. This 
effect will be considered in §4. 

Given that the thin line case has already been 
explicitly treated, we choose to focus on the thick 
line case. The major simplification that results 
comes in the escape probability. Recall that (3 
varies between and 1. For thin lines, (3 approaches 
the value of unity, meaning that photons scatter 
isotropically. For thick lines, j3 decreases to small 
values. Taking 77 > 1, the escape probability reduces 
to 



P 



2 

3T 



—J (l-^ 2 )^ (27) 
(28) 



Clearly at large radius where 77 drops below 
unity, the approximation breaks down. Defining 
7 = 3TR*/2r, Figure || shows a plot of 7 • (3, 
indicating that this product begins to plateau at 
unity around 7 of 2. Since the constant expansion 
approximation applies only for r/R* exceeding 
about 2-3, values of T J> 10 are required for lines 
to show interesting optical depth effects and for 
expression ( p8| ) to apply. 

Continuing in the assumption that the line 
is optically thick, and noting that for constant 
expansion e = €{)(R*/r) 2 , the source function 
reduces to the form 



S v (T)n\T€oB v ( — 

2 \ r 



(29) 



The emergent intensity from a given point (r, /1) is 
given by 



S v {r)e- T "^\ 



(30) 



where the continuum optical depth can be derived 
analytically (e.g., MacFarlane et al. 1991; Ignace 
2001) as 



r c (r,/x) = T 



sin# 



(31) 
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which is to be contrasted against the analogous 
asymptotic limit for thin lines from Ignace (2001): 



for tq(E) = k c (E)M /AttR^Voo, the total continuum 
optical depth to the star along the line-of-sight. 

To calculate the emission profile then requires an 
integration over the isovelocity zone corresponding 
to each frequency in the profile. For constant 
expansion, the isovelocity zones take the form of 
cones with v z = — v^^l, so that constant v z implies 
fixed ji. The integral to be evaluated is then 



F l {v z ) = 



2tt 

d 2 



-M r sin 2 9 dr, 



(33) 



where we have substituted for pdp = sin 2 9rdr, and 
^min equals i?* for v z < 0, but i?*/ sin# for v z > 
owing to stellar occultation. An important point 
here is that the integral is over radius, so that the 
angular dependence is a constant of the integration. 



With a change of variable u 
reduces to 



R*/r, the integral 



Fi(v z ) 



2TrRie B u f 3T\ . 2 

s sin ' 

D 2 V 2 



- Toue / sinB du. 



3nRl e B v 
D 2 



sin 



1 - e" 



(34) 



sh (8i5) 



where 9 = cos~ 1 (—v z /v 00 ), and u max = 1 for 
v z < and u max = sin 9 for v z > 0. Figure || 
shows normalized emission profiles at various 
tq. The case of tq = 1 indicates that the wind 
attenuation of X-rays is only mild, which invalidates 
our approximation that that X-rays emerge only 
where v r = v^; however, the tq = 1 case is 
included to demonstrate that the profile shape varies 
qualitatively little for a broad range of tq values. 

At large tq, the line profile reduces to one of 
fixed shape described by 



Fi(v z ) 



3irR 2 e B v 
D 2 



sin d 



(36) 



Fi(v z ) 



D 2 Vno 



sin ( 



(37) 



The profile morphologies of the thin and thick 
asymptotic forms are compared in Figure ^. Note 
that both are asymmetric, with peak emission 
occuring at blueshifted velocities. However, the thick 
line is much more nearly symmetric in appearance, 
and its peak emission occurs not far from line center. 
In fact for the asymptotic limit of Fi oc sin 3 9/9, the 
location of the peak and the FWHM value of the 
profile can be easily computed. Peak emission occurs 
at v z = -0.24 x Voo, and the FWHM is 1.26 x 
(alternatively, the HWHM is 0.63 x Voo), which 
holds for all thick lines with T ^> a few. However, 
although resonance line scattering does conserve 
photon number, photons could be destroyed by 
other processes since there is a longer "dwell" time 
for photons than would be the case for optically 
thin lines that scatter only once. The next section 
considers the influence of reabsorption effects for the 
net line emission and profile shape. 



4. OPACITY EFFECTS ON 
ABUNDANCE DETERMINATIONS 

Optically thick resonance zones not only alter 
the profile shape, they also modify the fraction 
of X-ray line photons that are absorbed by the 
background continuum opacity prior to escape 
from the wind. They do this in two distinct ways. 
First, the preference for azimuthal emission slightly 
alters the average column depth through which 
the escaping flux must pass, and this changes the 
equivalent width of the line profiles. Second, if the 
Sobolev length ^ so b is not negligible, as would be 
the case for lines that are highly broadened due 
to X-ray temperatures or microturbulence, then 
optically thick resonance zones can "trap" photons 
long enough to increase the total pathlcngth over 
which continuum absorption can occur (Hillicr 
1983). Both of these processes were calculated 
accurately by Hummer & Rybicki (1985). If these 
effects increase the reabsorption, then observed line 
equivalent widths will be reduced, and anomalously 
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low abundances might be inferred. Indeed, this is 
the interpretation given by Kitamoto et al. (2000) 
to the substantial shortfall in metal abundance 
inferred from optically thin fits to ASCA X-ray line 
observations of several OB stars. However, they 
did not apply the theory of Hummer & Rybicki 
(1985) in making their estimates, and here we show 
that they overestimated the degree to which photon 
destruction between scatterings in the line can 
reduce the observed line flux. 



4.1. The Equivalant Width for Optically 
Thick Line Emission 

To quantify the change in absorption due to 
the first effect above, note that optically thick 
lines in constant-expansion winds emit photons 
preferentially in directions azimuthal to the flow. 
Thus optically thick line photons emitted in 
generally outward directions (i.e., emitted toward 
the observer from the near side of the wind) will 
have an overall reduced chance of escape, since more 
of them are emitted from high impact parameters 
which experience greater continuum optical depth 
for a given radius. Conversely, the photons observed 
from the far hemisphere will be augmented on 
average, because then emitting more photons at 
higher impact parameters implies less occultation 
by the underlying wind. To ascertain which of 
these competing influences dominates, we write the 
expression for the ratio of the photon escape rate in 
the optically thick and thin limits, which also gives 
the equivalent-width (EW) ratio 



EWi 



thick 



/~ dr r 2 p 2 ft d9 1.5 sin 3 9 e -nR.e/rOn, 



EW X 



ihii 



f™dr r 2 p 2 ftd9 e -ToR*e/ri 



(38) 



In the above expression, equation ( |32| ) was used 
for the photoabsorptive optical depth, and the 
normalized emission profile from an optically 
thick resonance zone at constant expansion speed 
is 1.5 sin 2 9. It was also assumed that the line 
emissivity scales as p 2 , and that all line photons 
eventually escape their natal resonance region. Later 
we will consider the fraction that may be absorbed 
within the resonance zone. 



Applying the condition p oc r~ 2 for constant 
expansion allows the elementary radial integrals in 
equation ( ^8|) to be carried out explicitly, giving in 
the limit of large total radial wind optical depth tq 



EWi 



thick 
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= 1.5 
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(39) 



When this expression is evaluated numerically, 
it is found that optically thick angular emission 
actually enhances escape, so the far-hemisphere 
effect dominates. This is the opposite sense as would 
be needed to help explain the abundance anomalies, 
but is of little significance because the magnitude of 
the increase in the equivalent width is only 2%. 



4.2. Reduction in Line Equivalent Width 
By Continuum Reabsorption 

The second effect mentioned above, reabsorption 
within the finite region in which the X-ray line 
photons are created and scattered, will certainly 
reduce the equivalent width when the lines become 
optically thick. This destruction mechanism is 
assumed to be dominated by continuum absorption, 
rather than thermalization within the line itself, 
and was considered by Kitamoto et al. (2000) to 
be a promising solution to low equivalent-width 
anomalies in ASCA observations of several OB 
winds. Unfortunately, they assume that photon 
escape from the X-ray-line resonance zone is 
diffusive, which produces an overly large total 
pathlength during which reabsorption can occur. 
Instead, escape of photons from the line-scattering 
region is known to occur by a "single-flight" escape 
process (Hummer & Rybicki 1982), whereby Dopplcr 
redistribution enables photons to rapidly reach 
the line wings. At that point the mean- free-path 
approaches the scale of the region and the escape is 
non-diffusive. In such a process, the total pathlength 
traversed scales linearly with the length scale for 
escape, not quadratically as assumed by Kitamoto 
et al. (2000). 

Nevertheless, the total pathlength is indeed 
enhanced by optically thick line scattering, so some 
of the abundance anomaly may be explained by 
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this process. The basic theory of how continuum 
opacity affects line formation in the Sobolev 
approximation was developed by Hummer & Rybicki 
(1985), so only a rough sketch of their derivation 
is included here, expressed in intuitively motivated 
probabilistic terms. For a two-level atom in complete 
frequency redistribution, common assumptions for 
line formation in the Sobolev approximation, the 
line source function may be written 



eB 



- - + G ' (40) 

where e is the destruction probability per scattering 
and is assumed to be small, B is the Planck function, 
and £ gives the fraction of radiative downward 
transitions in the atom that are not balanced locally 
by radiative upward transitions. Since eB here 
describes the photon creation rate, and it is likely 
that e <C £, the essential statement for our purposes 
is that Si oc £ , due to multiple scattering within 
the optically thick resonance layer. Hence £ _1 
may be associated with the expected number of 
scatterings each photon undergoes between creation 
and escape. 

In the first-order escape probability 
approximation, which is especially applicable 
in the spirit of the Sobolev approximation, £ is 
replaced by the probability that a line photon which 
has just been emitted will not be scattered again in 
the line. If the only other possibility is escape, then 
in the notation above this implies £ = /3. However, 
the presence of continuum opacity allows for another 
alternative to scattering, which is absorption within 
the resonance zone. Thus £ is increased, and Si 
reduced, by the presence of continuum absorption 
(e.g., Puis & Hummer 1988). 

The expression for the probability that no line 
scattering will occur is most easily found by first 
writing the probability that a line scattering will 
occur. This is pieced together by first noting that 
the photon is assumed to be emitted isotropically, 
so is emitted randomly in a flat distribution over 
direction cosine fi. In complete redistribution, it will 
also be emitted randomly over the frequency profile 
4>{x), typically a Doppler profile. Given an emission 
within d/i of angle /i and within dx of frequency x, 
in the Sobolev approximation in an expanding wind, 
the photon propagates effectively not in real space 



but rather toward redder comoving frequencies x' . 

If the photon makes it to some x 1 along direction 
the probability of scattering within dx' of x' is 
given simply by T^4>(x')dx' , where r M is the Sobolev 
optical depth of the resonance zone along direction ^. 
In the absence of continuum opacity, the probability 
that the photon will make it to x' is thus given by 



P(x,x',fj,) = [*(*)-*(*') 



(41) 



where 



$(x) - ®{x') 



dy Hv) , (42) 



and it is canonical to set $(— oo) = via the 
definition 

<D(x) = f dy cf>(y) . (43) 

J — oo 

In all, one obtains 

1 f 1 f°° f x 
£ = /3=l-- / dfj, dx 4>(x) / dx t^P(x,x , 



(44) 

Here r M and (3 are given in equations ( f22[) and ( ^7| ) 
respectively, and if we denote the azimuthal Sobolev 
line optical depth as 

it* 



n = T- 

T 

then equation ( |44"| ) yields, in agreement with 
equation (E7D, 



3r, 



(45) 



(46) 



The introduction of continuum opacity simply 
reduces P(x,%' , [/,) by an additional factor which 
represents the probability of continuum absorption 
during propagation from x to x' . If the effective line 
optical depth over the interval dx' is dx l T ll 4>(x') 1 then 
the continuum optical depth over that propagation 



interval is dx'r^C, where 



C = 



Kl 



(47) 



and k c and ki are respectively the continuum 
and mean line opacities. By its definition, the 
continuum-absorption parameter C is readily shown 
to be related to the ratio of r c , the radial continuum 
optical depth to the resonance zone in question, 
and Ti, the azimuthal Sobolev optical depth of that 
resonance zone, via 
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in a steady spherical wind at constant speed 



Note once again that here Vth represents all of 
the velocity broadening over the Sobolev length, 
and so may include more than the thermal speed 
of the ion in question and may approach the 
hydrodynamical sound speed of the hot gas, which 
is at the level of tens of percent of the terminal 
speed. If Vth is much greater than about 1/10 of 
Uoo, the Sobolev approximation itself may come 
into question, so this may be viewed as a rough 
upper bound for the broadening to which our theory 
can apply. Nevertheless, since r c ~ 1 where the 
observable line emission forms, and r\ J> 10 is 
required to apply optically thick line theory, the 
condition Vth & Voo/10 implies that C <^ 10~ 2 . This 
unusually large upper bound for the continuum 
opacity, relative to the lines, is a manifestation of 
an implicit assumption that the hot line-forming gas 
represents only a small fraction of a wind dominated 
by cooler and continuously-absorbing gas. 

By integrating the continuum optical depth 
encountered by a photon propagating from x to x 1 , 
the absorption-modified expression for P[x,x',^i) 
then becomes 

P(x,x',fi,C) = e- Tl [*(*)-*(*')] x e -nc(x-x') ) 

(49) 

and so when £ is approximated by the new 
probability that a scattering does not occur, the 
result is 

A£(C) = 1 - \ ( dfx [ dx 4>(x) 

2 J-l J-oo 

x / dx'4>(x')T^P(x,x',fi,C) (50) 

J — oo 

This can be written 



where 



and 



A£(C) 



ac) = m + A£(c) (5i) 



^(0) = P = — (52) 
3tj 



1 poo 

d/i / dx <j)(x) 

J-oo 

x / dx'(j)(x')T fi P(x, x, n, C) (53) 



was derived in an alternate notation by Hummer & 
Rybicki (1985). 

The correction to £ is found to become 
significant when Cr; ~ 0.1, and when this is true, 
continuum reabsorption within the line approaches 
the static limit for an infinite atmosphere. Then 
the tabulated results of Hummer (1968) show that 
a reasonable approximation for A£(C) when C is 



within an order of magnitude of 10 is 
A£(C) = 4C . 



(54) 



This implies that the factor by which the line 
source function Si must be reduced to account for 
continuum absorption when Cti ~ 0.1 is 



(3 



1 



si(c) m_ 

5,(0) ' £(C) " (/3 + A£(C)) ' ' (l + 6Cn) ' 

(55) 

As mentioned above, Ct, ~ fth/woo near 
the X-ray photosphere, so when is of order 
10 3 km s _1 , Cn can reach values as large as 
0.1 only when the line broadening represents a 
substantial fraction of the sound speed in the X-ray 
gas. Interestingly, although the parameter C does 
affect the source-function correction at any given 
radius, it is not directly relevant to the integrated 
equivalent width, because increasing re c simply 
moves the resonance zones of interest out to larger 
radii and therefore lower density. This regulates 
the reabsorption within the resonance zone in a 
manner that is independent of C . Instead, it is the 
parameter uth/^oo which controls the correction to 
the equivalent width, so we give this parameter its 
own designation, u t h, and note that Cr, = Uth^c- 
The equivalent width correction is calculated 
by inserting Si(C)/Si(0) = (1 + from 
equation (55) into the integral in the numerator of 
equation (38), and the result is shown in Figure |^. 

It is apparent from this figure that the 
equivalent-width correction is only substantial 
when uth is appreciably elevated beyond the purely 
thermal motions of the heavy ions. This might 
occur, for example, if overpressures at the sonic scale 
in the hot gas were to accelerate flows stochastically 
on the scale of the Sobolev length, which would yield 
superthermal microturbulencc in the X-ray lines. In 
such a scenario, Figure || shows that the observed 
equivalent widths could be reduced by up to about 
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40% before the basic assumptions of this model 
are severely violated. Such a substantial reduction 
could certainly alter the inferred abundance of the 
ion involved, although perhaps not to the degree 
expected by Kitamoto et al. (2000). Whether 
or not it could provide a complete explanation 
for the observed abundance anomalies remains a 
question for more detailed physical models; the 
sole conclusion here is that if abundance anomalies 
are to be explained by continuum absorption 
within the line resonance zones, some type of 
microturbulent broadening is required to generate 
Sobolev lengths with Z so b ;> O.lr, straining the 
Sobolev approximation itself. This possibility is not 
beyond reason, however, as a radial variant on this 
same concept has already been invoked (Lucy 1982) 
to explain "black troughs" in UV resonance-line 
spectra of OB winds (Lamers & Morton 1976; 
Prinja, Barlow, & Howarth 1990). 

4.3. The Influence of Continuum 
Absorption on Line Profile Shape 

Appreciable values of C n not only effect the 
equivalent width in low-resolution spectra, they also 
alter the profile shape at high resolution. Since 
t c = To(R*/r) implies that Cn oc r^ 1 , the effect 
of continuum absorption on the line becomes less 
important at large radii. Thus when the factor given 
by equation ( |55| ) is inserted into equation (^3|), the 
profile shape is altered by the fact that Si will now 
fall less rapidly than the r~ 3 of equation (p9|). Since 
C is fairly small, this will only be important when 77 
is large, and in that limit the modified equation (|34| ) 
becomes 

(56) 

where we have defined 

r°° P -y 
f(x) = xe x dy . (57) 

Jx V 

This profile shape is plotted in Figure ^ for a 
range in values of wth/^oo- Note that if w t h *C «oo, / 
approaches unity and we recover the To 3> 1 limit of 
equation ([m|). The unphysical limit u t h 3> «oo can 



be used as a check, since then f(x) is of order x, and 
the profile approximates the (symmetric) parabolic 
shape 1 — v\ /v^ characteristic of /^-dependent line 
emissivity over a "skin depth" of fixed continuum 
mean-free-path. This limit is of course unphysical 
because it completely violates the approximations 
made, particularly the Sobolev approximation. 
Instead, the importance of the profile shapes seen 
in Figure ^ applies for much more modest values 
of Uth/Woo) which are seen to reduce the profile 
asymmetry even further than that due solely to 
continuum reabsorption outside the resonance zone. 

5. SUMMARY AND DISCUSSION 

Our paper considers the role of optical depth in 
X-ray emission-line formation in stellar winds. The 
line photons are assumed to be created by collisional 
processes, and thus scale with the square of the 
density, and may then be resonantly scattered. We 
have demonstrated that strong lines with sufficient 
opacity to multiply scatter emerging photons will 
alter the observed profile shape in a recognizably 
characteristic way. Our line synthesis has taken into 
account the "cool" wind photoabsorption, which 
can strongly attenuate X-ray photons at soft and 
intermediate energies. 

We argue that when the wind is highly optically 
thick to photoabsorption, it is adequate to consider 
constant expansion for the purpose of computing 
the line profile. In this limit the expressions for 
the radiative transfer can be solved analytically. 
Lines that are optically thick over a resonance 
zone are found to have profile shapes that scale 
as Fi (x sin 3 0/0, which is to be compared to 
the analogous thin- line case of sin 0/0, where the 
observed velocity shift is v z — —Voq cos 9. 

' It useful to compare these theoretical values 
against those measured for £ Pup by Cassinelli et al. 
(2001). These authors list the HWHM and ccntroid 
shift of five lines in their Table 2 (we ignore the 
N VII profile which is clearly different owing to its 
flat-topped morphology). The centroid shifts range 
between -460 and -700 km s~\ and the HWHMs 
range from 970 to 1340 kms -1 . Cassinelli et al. 
quote a terminal speed of 2250 km s _1 . Hence peak 
line emissions fall between —0.2 and —0.3 of Voo, and 
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the HWHM are around 0.4 to 0.6 of v x . This data is 
compared with the asymptotic thick line predictions 
in Figure [?] (i.e., without reabsorption effects). The 
HWHMs all fall below our predicted value, but the 
centroid values are all consistent with our thick line 
results. This is somewhat surprising in light of our 
rather simplified assumptions and the fact that we 
do not allow for any radial variation of the hot gas 
temperature. We must interject the caution that we 
have not actually tried to fit the observed profiles, 
only compare global characterizations in the form 
of where peak emission occurs and the overall line 
width. More detailed line fitting to the data is 
necessary to determine categorically whether thin 
or thick lines better represent each observation, but 
our simple comparison suggests that optical-depth 
effects could be a key element to interpreting the 
observations, and that the £ Pup data is more 
consistent with the standard wind-shock paradigm 
than previously thought. 

However for £ Ori, 9 1 Ori C, and S Ori, the 
observed emission profiles are shown or claimed to be 
extremely symmetric, with little or no centroid shift. 
(Note that for observed line profiles of modest S/N 
and spectral resolution, our theoretical line shapes 
could match the symmetric appearance only if the 
central line wavelengths were misplaced slightly to 
the blue.) Thick resonance-line scattering by itself 
is not able to explain such symmetries, at least not 
without allowing parameters such as Tx to vary with 
radius, or introducing new free parameters. However, 
a consideration of reabsorption by the continuum 
opacity within a highly optically thick resonance 
zone can in principle lead to symmetric lines with 
an inverted parabolic shape in the limiting case, and 
suppression of the line emission. The effect depends 
on the relative line broadening Wth/^ooj leading to a 
modest effect for only thermal broadening by the hot 
gas, but can significantly reduce the line equivalent 
width (« 40%) if there are extreme microturbulent 
motions of order 0.25voo. For this latter case, the 
influence of reabsorption by thick resonance-line 
scattering may help to understand both the line 
symmetry in these three stars and the sub-solar 
abundance problem. One major difficulty is whether 
observed lines have sufficiently high values of T so 
that our results apply. A second is, where are the 
optically thin lines that do show the theoretically 
expected morphology for distributed wind shocks, 



and do such lines also exhibit abundance problems? 

One effect that we have ignored has been 
the influence of the wind acceleration zone. In 
lower-density winds, the formation of lines at 
intermediate X-ray energies will sample these inner 
radial depths of the wind. This will lead to a source 
function of the form S v oc (p w /0) = T;/(r 2 u r ), where 
i> r varies rapidly in the span of just l-2i?* and 77 
depends on the velocity gradient. However, the wind 
attenuation will also be enhanced, plus reabsorption 
effects of line photons will be more severe. So there 
are strongly competing effects between shifting the 
profile to greater symmetry versus asymmetry, and 
a more detailed analysis will be required in general. 

Another effect that has been ignored is the 
influence of "porosity" . We have treated the hot 
X-ray emitting gas as being thoroughly mixed 
throughout the ambient wind flow. However, 
the wind is probably highly structured in density 
and velocity, both radially and laterally. If the 
porosity is severe, in the sense that cool absorbing 
wind material is strongly clumped such that the 
mean-free-path for X-ray photons becomes quite 
long, then this may explain the lack of profile 
symmetry in ( Ori, 6 1 Ori C, and 5 Ori (Cassinelli 
2001, private communication). Of course, one must 
then pose the question of why £ Pup is different. 

Finally, should optical-depth effects be relevant 
for other early-type stars? The B star winds are far 
less dense than O stars, so that except at the softest 
energies, the X-rays from the entire wind probably 
emerge free of attenuation. In this case, the line 
synthesis presented by Owocki & Cohen (2001) will 
likely be able to match the profiles, should they be 
observed. On the other hand, the Wolf-Rayet winds 
are substantially denser than for the O stars, so 
dense that n is 10's and even 100's of i?* in extent. 
In this case, strong lines are probably optically thick 
over much of the wind, but this is compensated by 
the fact that the observed line emission emerges 
from large radius. Such lines can certainly be 
assumed to form in a constant-expansion wind, so if 
thin, the asymptotic profile shape derived by Ignacc 
(2001) should apply, and if thick, the asymptotic 
profile shape derived here applies instead. These 
conclusions need to be tested by future observations. 
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Fig. 1. — A plot of the optical depth unity 
radius from the wind photoabsorptive opacity versus 
wavelength in the X-ray band. Typical O star and 
wind parameters are assumed. The horizontal dotted 
line indicates r\ = , and the vertical line gives the 
location of E = 1 keV, corresponding to A = Y1AA. 
The location of strong emission lines are noted. 




2 4 
7 = 3TR,/2r 



Fig. 2. — A plot of how rapidly the line escape 
probability approaches the asymptotic limit of (3 = 

7 _1 = 2r/3TR :t . The vertical is shown as the product 

This preprint was prepared with the A AS IATeX macros 7 P- Clearly (3 is nearing its asymptotic form at 7 « 2. 
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Fig. 3. — A comparison of normalized thick line 
profiles showing the influence of the wind attenuation 
with To = 1, 2, and 10 as indicated. As long as 
the line is optically thick in resonance line scattering, 
the width and centroid shift of the line peak are not 
strongly affected by the wind photoabsorption. 
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Fig. 4. — A comparison of thin versus thick X-ray 
emission lines, both in the asymptotic regime of large 
To- Both lines are asymmetric with a blueshiftcd 
centroid of the emission peak, but the thick line is 
more symmetric than is the thin case. 
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Fig. 5. — A plot of how the equivalent width (EW) 
of thick lines are reduced as a function of the degree 
of continuum reabsorption as parametrized by the 
line broadening u t h = wth/«oo- The normalization 
EW(0) corresponds to the line equivalent width with 
u t h = 0. In the presence of microturbulence, values of 
Uth may rise to a couple of tenths leading to a modest 
reduction (20-40%) in the line EW. 
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Fig. 6. — A plot of different thick line profiles as 
wth/foo is increased from to 0.16 (solid lines), with a 
thin line profile shown for comparison (dashed) . The 
influence of reabsorption is twofold: larger values of 
the broadening diminish the overall amount of line 
emission, and the profile becomes more symmetric 
in shape. In the unphysical limit of fth/foo 3> 1, 
the profile will take on the purely symmetric form of 

(l-«z 2 A>co 2 ). 
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Fig. 7. — A comparison of the HWHM and centroid 
shift (both normalized by Uqo) for £ Pup observed 
with Chandra versus the predictions of our limiting 
thick line case. The points correspond to five different 
lines with errorbars as tabulated in Cassinelli et al. 
(2001). Although the HWHMs are slightly smaller 
than predicted, the centroid values agree well with 
our theoretical results. Curiously, lines at longer 
wavelength appear more nearly consistent with our 
predicted result. At the longer wavelengths, the line 
emission escapes from larger radius (see Fig. |l|) where 
the constant expansion approximation that we have 
used is best satisfied. 



